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1. Introduction 

Refinable funetions are functions that are in a sense self-similar: If you add shrunken 
translates of a refinable function in a weighted way, then you obtain that refinable 
function again. For instance, see Figure 1 for how a quadratic B-spline can be de- 
composed into four small B-splines and how the so called Daubechies-2 generator 
function is decomposed into four small variants of itself. 

All B-splines with successive integral nodes are refinable, but there are many 
more refinable functions that did not have names before the rise of the theory of 
refinable functions. In fact we can derive a refinable function from the weights of 
the linear combination in the refinement under some conditions. 

Refinable functions were introduced in order to develop a theory of real wavelet 
functions that complements the discrete sub-band coding theory.^ Following the 
requirements of wavelet applications, existing literature on wavelets focuses on re- 
finable functions that are £2-iiitegrable and thus have a well-defined Fourier trans- 
form, are localized (finite variance) or even better of compact support. It is already 
known, that polynomial functions are refinable as well.^ In this paper we want to 
explore in detail the connection between polynomials and the respective weights for 
refinement. 

Our results can be summarized as follows: 



1 



31, 2012 3:6 WSPC/WS-IJWMIP paper 



2 Henning Thielemann 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Quadratic B-spline Daubechies-2 wavelet generator 

Figure 1. Refinement of a quadratie B-spHne and the orthogonal Daubechies-2 generator. The 
blaek line is the refinable function that is composed from shrunken, translated and weighted 
versions of itself, displayed with dashed lines. 

• Masks that sum up to a negative power of two refine polynomials that are 
uniquely defined up to constant factors. Other masks are not associated with 
a polynomial. (Theorem 2.1) 

• For every polynomial there are infinitely many refinement masks, and these 
refinement masks can be characterized in a simple form. (Theorem 2.2 and 
Theorem 2.3) 

• There is a simple iterative algorithm for approximating a polynomial that is 
associated with a mask. (Theorem 2.4) 

2. Main Work 
2.1. Basics 

We start with a precise definition of refinable functions. 

Definition 2.1 (Refinable function). The vector m with m G and a finite 
number of non-zero entries (m e io (^) ) is called a refinement mask for the func- 
tion (fi if 

^{t)=2-J2mj-V>{2-t-j) (2.1) 

holds. Vice versa the function if is called refinable with respect to the mask m. 

The factor 2 before the sum is chosen, such that the following law (Lemma 2.1) 
about convolutions holds. Unfortunately this enforces adding or subtracting 1 here 
and there in some of the other theorems. There seems to be no convention that 

asserts overall simplicity. 

Definition 2.2 (Convolution). For sequences h and g the convolution is defined 
by 

[h *g)k = ^ hj ■ gk-j 
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and for real functions the convolution is defined by 

{(fi*ip){t)= / (p{t) ■ ilj{t — t) dr 
Jm. 

Lemma 2.1. If (p is refinable with respect to h and ip is refinable with respect to g, 
then (p^tjj is refinable with respect to h* g. 

For a proof see 8. For the proof of our theorems we need two further lemmas 
about differentiation and integration. 

Lemma 2.2. // the function (p is refinable with respect to mask m, then its deriva- 
tive if' is refinable with respect to mask 2 • m. 



Proof. Derive both sides of the refinement equation with respect to t. 



ip{t) = 2-^mj-^{2-t-j) 
^'(i) = 2-^m,-2.(p'(2-f-i) 

Lemma 2.3. // function ip is refinable with respect to mask m and there is an 
antiderivative $ with = J^(f{T)dT, then the antiderivative i->- + c) is 
refinable with respect to mask | • m where the constant c must be chosen as follows: 

• For ^2- m,j ^ 1 it must be c = ^/^^ ''^ ■ 

• For ruj = 1 and TTij ■ = the constant c can be chosen arbitrar- 
ily. 

• For J2j f^j = 1 o.'i^d ruj ■ i) 7^ there is no valid value for the con- 
stant c. 

These choices for c are the only possible ones. 



Proof. We start with the necessary condition; that is, given that the antideriva- 
tive {t !->■ + c) is refinable with respect to mask 5 • m, what are the possible 
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integration constants? 

f2-t-i 



2 

iez 



$(t) + c-(l-^mj) • / (^(r)d^ 



^rnr(f ^(r)dT + $(-j) 



Y^ruj-U- I V'(2-r-j)dr + $(-j) 



= / |^E2.m, -y.(2-T-j)j d 
c ■ (1 - ^ m,) = ^ m,- ■ 



Proof of the sufficient condition: By substituting c by the admissible values we 
verify, that the antiderivative with that offset is actually refined by ^ • m. □ 

Remark 2.1. For ^jTrij = 1 and m and ip with support, that is bounded on 
at least one side, the second case of Lemma 2.3 applies. This means that all an- 
tiderivatives of if irrespective of the integration constant are refinable with respect 
to ^ • m. 

Sketch of the proof: Without loss of generaUty let m and ip have support, that 
is bounded at the left. The refinement equation implies, that the bounds of their 
support arc equal. If their support is entirely on the positive real axis, then Vt < 
^{t) = and further on in all summands of mj ■ j) at least one factor is 
zero. This implies mj ■ j) = 0. 

Now, when ip is refined by m and k is an integer, then (p translated by k is 
refinable with respect to m translated by k. This way we can reduce all m and ip 
with bounded support to ones with support on the positive real axis. 



2.2. Conversions between polynomials and masks 

If we generalize refinable functions to refinable distributions, then the DiRAC im- 
pulse is refined by the mask S with Sj = < and the k-th derivative of the 

' lo :i^0 

DiRAC impulse is refined by 2''-5. Vice versa the truncated power function (t i-^- f^) 
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(t :t>0 ^ ^ 
with G N and t+ — < is refined by 2 ^ ■ o. Intuitively said, truncated 

+ \0 :t<0 

power functions are antiderivatives of the Dirac impulse. 

Once we are thinking about truncated power functions, we find that ordinary 

power functions with natural exponents are also refinable. Then it is no longer a 
surprise, that polynomial functions are refinable, too. For example / with f{t) = 
1 + 2 • i + is refined by the mask i ■ (3, -3, 1): 

2 . 1 . (3 . f{2t) - 3 • f{2t - 1) + 1 ■ f{2t - 2)) 

= ^ • (3 • (1 + 2 ■ 2t + (2t)2) - 3 . (1 + 2 ■ (2t - 1) + (2t - 1)^) 

+(1 + 2- (2i-2) + (2i-2)2)) 
= ^ • (3 • (1 + 4i + 4*2) - 3 • 4*2 + (1 - 4i + 4*2)) 
= 1 + 2 ■ * + i2 . 

Now that we have an example of a refinable polynomial function, we like to 
know how we can find a mask that refines a polynomial function. Vice versa we 
want to know a characterization of masks that refine polynomial functions and 
what polynomial functions can be refined by a given mask. 

Before we start answering these questions we would like to stress the difference 
between a polynomial and a polynomial function. 

Definition 2.3 (Polynomial and Polynomial function). A polynomial p of 
degree n is a vector from Ri*^'- ""} . We need this for the actual computations and for 
performing linear algebra. A polynomial function pis a real function. The refinement 
property is a property of real functions. The connection between polynomial and 
polynomial function is 

n 

p{t) = j2pk-t'' . 

Our first theorem answers the question, "What polynomial can be refined by a 
mask?" 

Theorem 2.1. Given a mask m that sums up to 2~"~^ for a given natural num- 
ber n, there is a polynomial p of degree n such that m refines p. With the additional 
condition of the leading coefficient being 1, this polynomial is uniquely determined. 

Proof. We show this theorem by induction over n. 
• Case n = 

We want to show that a mask m with sum i can only refine a constant 
polynomial. Thus we assume contrarily that m refines a polynomial with a 
degree d greater than zero. In the refinement relation 

p(t) =2.^m, •p(2.*-j) 

jez 



31, 2012 3:6 WSPC/WS-IJWMIP paper 



6 Henning Thielemann 

we only consider the leading coefficient, that is, the coefficient of t''. 

= 2<Pd 

From d > it foUows that Pd = 0- 

Thus the degree d must be zero and by normalization it must be po = 1- We 
can easily check that this constant polynomial is actually refined by any mask 
with sum i. 
• Case n > 0: Induction step 

The induction hypothesis is that for any mask with coefficient sum 2~" we 
can determine a refining polynomial of degree n — 1, that is unique when 
normalized so that the leading coefficient is 1. The induction claim is that 
for a mask m with sum 2~"~^ we have a uniquely determined polynomial of 
degree n with leading coefficient 1. We observe that 2 -to satisfies the premise of 
the induction hypothesis and thus there is a polynomial q of degree n— 1 that is 
refined by 2 ■ m and that is unique when normalized. Since the coefficient sum 
of 2 • rn is at most |, it is diflferent from 1 and thus the first case of Lemma 2.3 
applies. It lets us obtain the m-refinablc polynomial p in the following way: 
Let Q be the antiderivative polynomial of q where the constant term is zero, 
then it is 

Vfc > = ^ 

^°= Q„-(l-2-») -g"^^-^(-^'^ • 

Now we turn to the question of why p is uniquely determined. Assume we 
have two normalized polynomial functions po and pi that are both refined by 
mask m. Then their derivatives po' and pi' are refined by mask 2 • m. Due to 
the induction hypothesis the normalized polynomial functions of po' and pi' 
are equal. The first case of Lemma 2.3 implies that the antiderivatives with 
respect to po' and pi' have the same integration constant, and thus po = p\. 

□ 

Theorem 2.2. Given a polynomial p of degree n, there is a uniquely defined mask m 
of support in {0, . . . , n} that refines p. 

For the proof of that theorem wo introduce some matrices. 

Definition 2.4. We express shrinking a polynomial by factor k by the matrix Sk- 

Sk e m{0--"}x{o.-,"} 
5fc = diag(l,fc,...,fc") 
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We represent translation of a polynomial by 1 by the matrix T and translation of a 
distance i by the power T' . 

T e ]]j{0'---,n}x{0,...,n} 

[0 ■.3>k 
The proof of Theorem 2.2 follows. 



Proof. We define the matrix P that consists of translated polynomials as columns. 

Now computing m is just a matter of solving the simultaneous linear equations 

p = 2- S2Pm 

We only have to show that P is invertible. We demonstrate that by doing a kind 
of LU decomposition, that also yields an algorithm for actually computing m. Our 
goal is to transform P into triangular form by successive subtractions of adjacent 
columns. We define 



what satisfies 



Ap = Tp-p 



A{T^p) = T'^Ap 



In the first step we replace all but the first columns of P by differences, yielding 
the matrix Ui. 



L7' = 



Ui = 


{p,Ap,TAp,...,T^-'Ap) 














n -1 


••■ 0\ 




nil-. 


• ^\ 


1 


_1 ... 




1 1 •• 


■ 1 





1 ■•• 


Li = 


1 •• 


■ 1 


\0 


••■ l) 




\0 • • 


• 1/ 



In the second step we replace all but the first two columns of Ui by differences (of 
the contained differences), yielding the matrix 1/2- 

U2 = ip, Ap, A^p, TA^p, T^-^A^p) 



31, 2012 3:6 WSPC/WS-IJWMIP paper 



8 Henning Thielemann 



/I • 


• 0\ 


1-1- 


• 


00 1 


. 





/I •• 


• 0\ 




1 1 •• 


• 1 




1 ■• 


• 1 




\0 ■ • 


■ 1/ 



We repeat this procedure n times, until we get 

C/„ = (p,Ap,A2p,...,A» 

P = Un ■ Ln L2 • Li 

Since the fc-th difference of a polynomial of degree n is a polynomial of degree n — k, 
the matrix t/„ is triangular and invertible. Thus P is invertible. We get 

= \-L-^-L-^ L-^ ■ ■ S.p , 

where the product C/,^^ ' ^^P can be computed by back-substitution and the multi- 
plications with LJ^ mean computing some differences of adjacent vector elementO 

Theorem 2.3. If p is a polynomial of degree n and in is a mask that refines p, 

then for every mask v the mask m-|-w*(l,— 1)"+^ refines p as well. Only masks 
of this kind refine p. The expression (1,— l)""*"^ denotes the {n + l)-th convolution 
of the mask 5o - 5i, that is (1, -1)° = (1), (1, -1)^ = (1, -1), (1, -1)^ = (1, -2, 1) 
and so on. 



Proof. We denote the convolution of a mask m with a polynomial by the ma- 
trix Cm- 
Cm e R{0.--"}x{0,...,n} 

Cm = Y,mi-T' 
The refinement equation can be written 

P = 2 • S2Cm ■ P 

Since Cm is a Laurent matrix polynomial expression with respect to T, it holds 

(^h+g = Ch + Cg 
Ch*g = Ch ■ Cg 
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2 • S'2Cto+„*(1,-1)" + i • p 

= S2{2-C^-p) + 52(2 • • p) 

= p + 52(2-a-(C("+_\)-p)) 

(n + l)-th difference of an n-degree polynomial vanishes: C^^]^^^ ■ p = 
= P 

We still have to show, that refining masks of p always have the form m + 
V * (1,-1)"+"^. Consider a mask hi that rcfiucs p. By computing the Laurent 
polynomial division remainder with respect to the divisor (1,-1)"+^ we can reduce 
hi to a mask ho that has support in {0, . . . , n} and we can reduce m to a mask g 
with support in {0, ... , n}, too. 

m = g + vo*{l,-lT+^ 
hi = ho + VI* {1,-1)"+^ 

Prom the above considerations we conclude that both g and ho refine p and the 
uniqueness property in Theorem 2.2 eventually gives us g = ho, thus 

hi=m+{vi-vo)*il,-lT+^ . □ 

Remark 2.2. By adding terms of the form v * (1, —1)"+^ we can shift the support 
of a mask and still refine the same polynomial. 



2.3. Example 

Por better comprehension of the theorems of the previous section let us examine a 
longer example. We would like to illustrate that the same refinement mask can refine 
different functions. However, if we restrict the function class to, say, continuous 
compactly supported functions or to polynomial functions, then a refinement mask is 
associated with a unique function. We would like to compare a continuous compactly 
supported function with a polynomial function, both being refinable with respect 
to the same mask. Unfortunately this is not possible, since for the former type of 
functions we need masks with sum 1, whereas for polynomial functions we need 
masks with sums that are powers of two that are smaller than 1. So, we are going 
to compare antiderivatives of the quadratic B-spline and polynomial functions that 
are refinable with respect to the masks jq ■ (1, 3, 3, 1), ^ • (1, 3, 3, 1), • (1, 3, 3, 1) 
according to Lemma 2.3. 

Pigure 2 shows the refinable functions. Since the quadratic B-spline as in Figure 1 
is refinable with respect to the mask ^ • (1, 3, 3, 1), its antiderivative as in the top- left 
plot in Pigure 2 is refined by ^^-(1,3,3,1). The sum of this mask is ^ and thus the 
only refinable polynomial for that mask is a constant polynomial. We normalize it 
to 1. This is shown in the top-right plot of Pigure 2. 
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Figure 2. Refinement of antiderivatives of the quadratic B-spline (left column) and polynomial 
functions (right column) with respect to the masks jg ■ (1,3,3,1),^ • (1,3,3,1),^ • (1,3,3,1) 
(first to last row). The meaning of the lines is the same as in Figure 1. 



The second row is associated with mask ■ (1, 3, 3, 1). We get the refined poly- 
nomial by integrating the function t 1 as in the proof of Theorem 2.1: 



Q{t) = t 

Pi = 1 = • ^ • (1 • + 3 • (-1) + 3 • (-2) + 1 • (-3)) = -\ 

3 



That is, t H^- — I + Hs refined by ^ ■ (1, 3, 3, 1). We repeat this procedure in order 
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to get to the last row of Figure 2. We start with the antiderivative of 1 — | + t: 

P. = l n = -S Po=^y^-^-(l-0 + 3.2 + 3.5 + 1.9) = ^ 
= ^ - 3 • i + 

We want to check whether the obtained polynomial function t |— 3-t + f^ 
is actually refined by ^ • (1, 3, 3, 1) according to (2.1): 



p{2t- 


0) = 


5 

2 " 




1 


p{2t- 


0) = 


5 

2 " 


-e-t + A-t^ 


p{2t- 


1) = 


13 
T 


-IQ-t + A-t^ 


3 


p{2t- 


1) = 


39 
Y 


-SO-t + 12-t^ 


p{2t- 


2) = 


25 
T 


-lA-t + A-t^ 


3 


p{2t- 


2) = 


75 
Y 


- 42 • i + 12 • 


p{2t- 


3) = 


41 

T 


-l^-t + A-t^ 


1 


p{2t- 


3) = 


41 

Y 


- 18 • t + 4 • 
















80 


-96-t + 32-t^ 



^■{80-96-t + 32-t^) = ^-3-t + f 
64 2 

That is, the mask actually refines the polynomial function. 

In the next step of our example we want to determine refinement masks for our 
polynomial functions by means of Theorem 2.2. However, it is already clear that we 
will not get the mask ^-(1,3,3,1) as a result because Theorem 2.2 only promises 
a mask of size 3 for a quadratic polynomial. Nonetheless this smaller mask should 
be compatible to the original one in the sense of Theorem 2.3. 

= ^ - 3 • i + 

Ap{t)=p{t-l)-p{t) 

= g-3.(t-l) + (t-l)2)-Q-3.t + t^) 

= 4-2-i 
A2^(t) = Ap{t - 1) - Ap{t) 

= (4-2.(t-l))-(4-2.i) 
= 2 
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I 4 



U2 = [p, Ap, A'p) =1-3-20 
1 Oy 

m=l . L^' ■ L^' - U^' ■ S^p = L^' . L^' - U^' - l ■ 



'10' 

-6 



32 



32 



1 



We leave it to the reader to verify that this mask actually refines our polynomial. 

The last thing we want to check is, that the difference between the short mask 
and the original mask is a convolutional multiple of (1,-1)^ as stated by Theo- 
rem 2.3: 

1 . (1, 3, 3, 1) - 1 ■ (1, 0, 3, 0) = 1 ■ (-1, 3, -3, 1) = -1 . (1, -If . 



2.4. The cascade algorithm 

We want to close the section on the theoretical results with an alternative way 
to compute the polynomial that is refined by a mask. It is an iterative algorithm 
for the approximate computation of a polynomial, and it is analogous to the cas- 
cade algorithm known for refinable functions of bounded support.^ The refinement 
relation 

p = 2-S2Cm-p 
is interpreted as a recursively defined function sequence with 

Pj + l = 2 ■ S2Cm ■ Pj ■ 

This iteration is in fact the vector iteration method for computing the eigenvector 
that corresponds to the largest eigenvalue. 

Theorem 2.4. Given a mask m that sums up to 2~"~-^ for a given natural num- 
ber n and a starting polynomial po of degree n that is not orthogonal to the refined 
polynomial, the recursion 

Pj+l = 2 • S2Cm ■ Pj 

converges and the limit polynomial limj_>.(x, Pj is refined by m. 
An appropriate choice for po is {0,. . . , 0, 1). 
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Proof. The matrix S^Cm expands to 




rrii : j <k 
:j>k 



and thus is of upper triangular shape. This implies that the diagonal elements 2^ ■ 
^i-i, TTii for j € {0, . . . ,n} are the eigenvalues. Because the mask sums up to 2~'^~^ 
the eigenvalues of 2 ■ S2Cm are {1, . . . , 2~"}. That is, the largest eigenvalue is 1 and 
it is isolated. These are the conditions for the vector iteration method, consequently 
the iteration converges to a vector that is a fixed point of the refinement operatioffl 

Remark 2.3. The eigenvectors of the eigenvalues are the respective derivatives of 
the main refinable polynomial. 

3. Implementation 

The presented conversions from masks to polynomials and back are im- 
plemented in the functional programming language Haskell as module 
MathObj . Ref inementMask2 of the NumericPrelude project^. However, note that 
the definition of the Haskell functions slightly differ from this paper, since the fac- 
tor 2 must be part of the mask. 

4. Related work 

So far, refinable functions were mostly explored in the context of wavelet theory. 
In this context an important problem was to design refinement masks that lead to 
smooth finitely supported refinable functions. It was shown that smoothness 
can be estimated the following way: Decompose the refinement mask m into the 
form (1,1)" *v, where n is chosen maximally. According to Lemma 2.1 this corre- 
sponds to a convolution of functions. However, strictly speaking, v corresponds to 
a distribution. The exponent n represents the order of a B-spline and is responsible 
for the smoothness of the refinable function, whereas for v there is an eigenvalue 
problem, where the largest eigenvalue determines how much the smoothness of the 
B-spline is reduced. 

The cascade algorithm^ was developed in order to compute numerical ap- 
proximations to refinable functions. A combination of the cascade algorithm and 
Lemma 2.1 was used by 1 for computing scalar products and other integrals of prod- 
ucts of refinable functions. This is required for solving partial differential equations 
using a wavelet Fajiepkhh approach. 

Discrete wavelet functions in a multiresolution analysis are defined in terms of 
refinable functions, but were not considered refinable functions at first. However in 
7 it is shown, that wavelets are refinable with respect to infinite masks. The trick 
is to use polynomial division for dividing the wavelet masks of adjacent scales: If y 
is refinable with respect to mask h, and V' is a wavelet with respect to mask g and 
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generator ip, then -ip is refinable with respect to ^ • h, where g t 2 is g upsampled 
by a factor of 2. 

Although polynomial functions are not in the main focus of the research on 
refinable functions, we got to know several discoveries of this relation after uploading 
our work to arXiv. The thesis 4 is the first reference known to us that explains how 
to obtain a polynomial function that is refined by a mask. The approach in this 
thesis is also based on polynomial differentiation. 

The first source known to us that describes the opposite way, i.e. how to find 
a refining mask for a polynomial function, is 3.^ In terms of our matrices the au- 
thor of that article does not just use the matrix P of shifted polynomials, but its 
factorization P = A-V. The matrix F is a Vandermonde matrix with 

V G ]R{0.---."}x{Ov,"} 

The matrix A is defined using the derivatives of p by A = , or 

element-wise by 

A G ]R{0,---,n}x{0,...,n} 

\{'V)-Pi+3 ■i + j<n 
I :i + j>n 



The Taylor expansion of p allows to express translations of p in terms of its 
derivatives: 

Thus multiplying A and V yields the sequence of shifted polynomials with x G 
{0,-1,..., -n}. 

With this factorization the invertibility of P follows obviously from the invertibil- 
ity of the triangular matrix A and the invertibility of the Vandermonde matrix V 
with respect to pairwise distinct nodes. 

Actually, in the article the general VANDERMONDE matrix with pairwise distinct 
nodes {£o, ■ ■ ■ ,^n} is used: 



V = 



/I 1 •■• 1 \ 



V(-£o)" Ml)" ••• {-^uTJ 



That is, a refinement mask for a polynomial can also be found when only a certain 
set of n + 1 nodes is allowed to be non-zero. The nodes may even be non-integral. 

''It is written in 2011, but was already presented in talks in 2003. 
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In retrospect we could have conducted our proof of Theorem 2.2 with arbitrary 
nodes, too. We would have to define P = {T^° ■ p,. . . , T^" ■ p) and then use divided 
differences instead of simple differences for the LU decomposition. 

The paper 2 extends the previous one by an exploration of the refinability of 
rational functions. They find that a rational function ip is refinable if and only if 
there is a real sequence s (s G (Z), i-e. a Laurent polynomial) and a positive 
natural number k such that 

that is, 3fh{z) G m{z) ■ 's{z) = 5(2;^), where 2^°"^ • m is the refinement 

mask. 

5. Future work 

There arc some obvious generalizations to be explored: refinement with respect to 
factors different from 2, separable multidimensional refinement and most general 
multidimensional refinement with respect to arbitrary dilation matrices. 

Another interesting question is the following one: By Lemma 2.1 we know, that 
convolution of functions maps to convolution of their refinement masks. We can use 
this for defining a kind of convolution. In order to convolve two functions (fo and (fi , 
we compute refining masks mo and mi , respectively, convolve the masks and then 
find a function that is refined by mo *mi. In case of polynomial functions there is no 
notion of convolution because the involved integrals diverge. We can however define 
a convolution based on refinement. Unfortunately, the mapping from a polynomial 
function to a refinement mask is not unique, consequently the defined convolution 
is not unique as well not to speak of the arbitrary constant factor. If we choose; 
arbitrary masks from the admissible ones, then the convolution is not distributive 
with addition, i.e. tp * {(po + <fi) = ip^ipo + ^pfipi is not generally satisfied. The open 
question is, whether it is possible to choose masks for polynomials, such that the 
polynomial convolution via refinement is commutative, associative and distributive. 
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